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Abstract: We demonstrate a holographic reahzation of color-flavor locking phase, 
using A/" = 4 SU{Nc) SYM coupled to Af = 2 Nf fundamental hypermultiplets as an 
example. The gravity dual consists of Nc D3-branes and Nf D7-branes with world 
volume gauge fleld representing the baryon density. Treating a small number Nc C 
Nc of D3-branes as Yang-Mills instantons on the D7-branes, we consider possible 
potential(s) on their moduli space or equivalently the Higgs branch. We show that a 
non-trivial potential can be generated by including the backreaction of the baryonic 
density on the D7-branes, this dynamically drives the instantons (= D3-branes) into 
dissolution. We interpret this as a color-flavor locking since the size of the instanton 
is the squark vev, and study the symmetry breaking patterns. Extending to flnite 
temperature setup, we demonstrate that color-flavor locking persists, and the thermal 
effect provides additional structures in the phase diagram. 



Contents 



i 



Introduction and Summary 
Instanton on the Flavor Branes 



2.1 



2.2 



Review of tlie D3D7 System at Finite Baryon Density 
Instanton Solution 



27^ Conformal Metric and Explicit Instanton Configuration 



1^. Linearized Supergravity Backreaction 



3Tl| Smeared Baryon Vertices 

3^ Backreaction from the D7-brane Electric Flux 



Dissolution of the Instanton and Color-Flavor Locking 

0| Additional Chern-Simons Term 



4^ Potential on the Instanton Moduli Space 
O What is Locked? 



|5|. Extension to Finite Temperature System 



STT] The D3D7 System at Finite Baryon Density and Temperature 
5^ CS Term from Backreaction 



|6|. Discussions 

1^. Check of Consistency for the Linearized Perturbation 



i 
i 

I 

n 

12 

El 
u 

17 

m 
n 

21 

23 
^ 
27 

29 

|3T| 



1. Introduction and Summary 

In QCD phase diagram, the color-flavor locking (CFL) phase, or more generically, 
the color superconducting phase, is expected to be present in a region with large 
chemical potential fi for baryon number. Perturbative analytic study of this phase 
(see Ref. for reviews) has mainly been done for very large fi such that the QCD 
coupling is weak. However, the issue on possible phase transitions from the hadronic 
phase at finite fi has not been addressed, as the system becomes strongly coupled. 
So far, neither direct experimental search, nor the lattice QCD simulation with "sign 
problem" have reached such region in the phase space. 
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Holographic techniques from gauge/string duahty 0, ^ may offer new insights 
to such issue, as they enable us to probe the strongly coupled region(s) in the phase 
diagram for QCD-like theories. Although the duality strictly works for large number 
of colors Nc':^ 1, the holographic techniques applied to QCD-like theories (so-called 
"Holographic QCD" ) have been rather successful in reproducing qualitative and semi- 
quantitative features of low energy QCD dynamics. In this paper, among other 
things, we shall show that a color-flavor locking occurs for a toy QCD-like theory at 
zero temperature, when the baryon chemical potential exceeds its critical value. ^ 

Problems of CFL in Holographic QCD 

Before entering the details on how to realize CFL phase in our model, let us 
summarize here the possible difficulties in obtaining it in holographic QCD. 

• In gauge/string duality, to treat Nf flavor branes as probes, we typically need 
to take Nc ^ Nf, while the CFL refers to a locking of the SU{3) flavor and 
the SU{3) color symmetries, i.e. Nc = Nf. 

• In gauge/string duality, usually only gauge-invariant quantities are considered, 
while in the CFL phase the order parameter is a gauge variant di-quark con- 
densate. 

The flrst problem is strictly technical, as when Nf ~ Nc, the backreaction of the 
flavor branes cannot be ignored, and renders it difficult to analyze in supergravity.^ 
Our approach used in this paper is to flrst separate some flnite number of color branes 
Nc{i.e. Nc ^ Nc), and investigate the locking of SU{Nc) color symmetry with the 
flavor symmetry. Though this procedure of separation is artiflcial, our result may 
suggest a piece of the whole picture. Another concern for the flrst problem is that 
in the strict Nc ^ oo limit, the theory does not reveal the CFL phase [§]. We don't 
consider this concern, since we will not perform a comparison with the chiral density 
wave (which is supposed to be favored at the strict Nc ^ oo limit) in our toy model, 
and also because a large but flnite value of Nc may give the CFL phase even for the 
real QCD 0. 

As for the second problem above, it is familiar to us that gauge-invariant corre- 
lators of QCD-like theories can be computed in their gravity duals, but in fact there 
are some gauge-covariant quantities which one can also compute in the gauge/string 

^Disclaimer: Note that our theory is not QCD but rather a supersymmetric generaUzation of it, 
and we shall only treat the squark condensation for the CFL. For a field-theoretical treatment of 
the squark condensation, see for example Rcf. Q]. 

^There are examples in which fully backreacted geometry is obtained (see for example Ref. 
and it would be interesting to generalize our results to those examples. 
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Figure 1: The structure of the phase diagram of the M 
is not accurate. See |l^, 21] for details.) 



2 SQCD. (Scales in this figure 



duality."^ In this paper, we use holographic techniques for Coulomb phase of super- 
symmetric Yang-Mills (SYM) theories [jlOl , pjj] , where a part of the gauge symmetry 
decouples from the rest. When the rank of this decoupled gauge subgroup is small, 
we may treat them in the same way as the probe flavor branes, and their gauge 
symmetry is manifest in the dual gravity description. We shall describe this in detail 
later. 

Supersymmetric QCD, the Holographic Dual and Phase Diagram 

The toy model we shall focus on is A/" = 4 SYM coupled to A/" = 2 fundamental 
matter hyper multiplets. The holographic dual of this theory was proposed by Karch 
and Katz |]12|, as a minimal deformation of the Af = 4 SYM to include quarks. 
The quark superfields are introduced as the lowest excitation on a string connecting 
Nc D3-(color-)branes and Nf D7-(flavor-)branes. For Nc ^ Nf, D3-branes can be 
replaced by AdS^ x geometry, and the flavor dynamics of the strongly coupled 
large Nc SQCD can be analyzed by the probe flavor D7-branes in that geometry. 
The quark mass m quark is proportional to the distance between the D3-branes and 
the D7-branes. 

For zero temperature T = and /i = 0, quarks and gluons are deconfined while 
quarks can form deeply bound mesons. The phase structure of this theory has been 
analysed by the holographic duality [1^ |T9|, ^ and at the leading large 
Nc expansion, it is known that there are two phases in the (/i, T) diagram: the 
meson phase and the melted meson phase (see Fig. |ip. In both phases, gluons are 
deconfined. In the meson phase, quarks are bound to form mesons with their discrete 
spectrum,^ while in the melted meson phase, the meson spectrum is continuous, and 



^Examples of that kind include computations explicitly uses string worldsheets in the dual 
gravity backgrounds; gluon scattering amplitudes, drag forces, quark-antiquark forces, heavy meson 
spectroscopy and Regge trajectory. 

^The meson spectrum at zero baryon density is studied in Refs. 
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Figure 2: Two embeddings of the D7-branes in the geometry. The shaded ball denotes a 
black hole with a horizon of the topology S^. Left: Minkowski embedding (meson phase). 
Right: black hole embedding (melted meson phase). 



there appears nonzero baryon number density. These two phases are characterized 
by the shape of the probe Nf D7-branes |9|, |2|, |6|, l^, m. For the finite 



temperature, the background geometry is an AdS black hole. The meson phase 
corresponds to the D7-branes away from the horizon, which is called "Minkowski 
embedding". On the other hand, in the melted meson phase, the D7-branes touch 
the horizon (see Fig. |^), and is called "black hole embedding". 

Since the local gauge symmetry U{Nf) on the D7-brane is identified as the global 
U{Nf)Y symmetry of the SQCD via the gauge/string duality, the chemical potential 
fi is identified as the value of the temporal component of the overall U{1) gauge field 
on the coincident D7-branes. In the meson phase, this gauge field is just a constant 
/X, while in the melted meson phase, there appears electric flux on the D7-branes: this 
configuration has a lower free energy which the holographic dual can compute, and 
thus favored. The phase transition is first order, and the critical chemical potential 
for T = is ficr = fn- In the melted meson phase, the shape of the D7-branes is 
a spike whose tip is inside the horizon. Electrically charged spikes on D-branes are 
identified as fundamental strings, so the existence of the electric flux means that the 
quark number density is nonzero. These are briefly reviewed in Sec. 

Dynamically Driven CFL 

Let us explain how the CFL phase of this theory can be realized in its gravity 
dual. First of all, note that our theory is A/" = 2 SQCD which includes squarks 
carrying the baryon (quark) number. So, once the chemical potential becomes large 
enough, we expect squark condensation, instead of di-quark condensation. We shall 
see this squark condensation in this paper: this is certainly a CFL, but also a Higgs 
phase. 
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As suggested before, we separate Nc D3-branes among N^. and treat them as 
probes, Nc <^ Nc- The relevant quark/squarks are strings connecting the Nc D3- 
branes and the Nf D7-branes. Condensation of strings connecting Dp-branes and 
D(p + 4)-branes is well-known [Q: the Dp-branes are dissolved into the D(p + 4)- 



branes, and the Dp-branes can be seen as finite size instantons on the D(p+4)-branes. 
Therefore, the CFL Higgs phase of the SQCD is equivalent, via the gauge/string 
duality, to the situation where the size of the instantons on the probe D7-branes is 
driven to become larger. We will show in this paper that this is indeed the case, by 
computing the potential of the instanton size modulus on the D7-brane U{Nf) gauge 
theory, in the melted meson phase. D3-branes are moved onto the D7-branes and 
dissolve on the D7-branes dynamically. 

This Higgs phase for T = was described in Refs. [^, |32|, |3^ (see also Ref. p^), 



and the potential for the instanton size modulus was considered in the absence of 
baryon density. At T = 0, the resultant potential vanishes (we review it in Sec. P]^ 



and |2.3D , thus there is no CFL. Our new point is that including a backreaction 
from the D7-brane electric flux (Sec. this generates a nontrivial potential for the 
intanton size modulus (Sec. |). The new potential has a run-away behavior (see 
Fig. H), causing the instantons to expand, hence the CFL Higgs phase is prefered. 
This new potential exists only in the melted meson phase, so, for T = 0, above the 
critical baryon chemical potential, the CFL Higgs phase appears — this is what we 
show using the gauge/string duality for the SQCD. 

The way this new potential emerges is quite intriguing. This is essentially a 
Chern-Simons (CS) term on the D7-branes, J ti F A F A F A C2. The backreaction 
of the electric flux on the D7-branes generate a nonzero constant Ramond-Ramond 
(RR) 3- form flux F3 = dC2 (|3.1| ).^ Substituting this to the CS term, we obtain 
J ti A A F A F, thus the electric potential At on the D7-brane interacts with the 
instanton density tiF AF, which gives a nontrivial potential ( [4.9| ).^ This generation 
of -F3 can also be thought of being sourced by baryon vertices, which are nothing 
but D5-branes wrapping (see also Ref. ||^). If one smears them, they provide 

a constant magnetic flux F3 along x^-x'^-x^ directions (Sec. |3.1D . So, our work is an 
example of backreacting baryon vertices. 



We also analyze the thermalized case with T 7^ (Sec. ^). Ref. showed that 



for the finite temperature, a nontrivial potential ( |5.23| ) for the size modulus of the 
instanton on the D7-branes is generated, before including the baryon density. This 
potential is minimized at a finite value of the instanton size. Therefore a Higgs CFL 
phase is prefered. Introducing baryonic density, we analyze the backreaction and our 



^The importance of this coupling between the NSNS 2-form and the F3 for the baryons was 
found in Ref. [||. 

similar CS mechanism was used for treating baryons Q in Sakai-Sugimoto holographic 
model ijsj, but used in a rather different way: the CS term was to stabilize the size of a single 
baryon in the model. 
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Figure 3: The structure of the phase diagram of the N = 2 SQCD given by a naive 
addition of the instanton size potential coming from the back reaction of the geometry. 
The melted meson phase (corresponding to the black hole embedding of the D7-branes) is 
devided into two distinct phases. The upper half denoted as "thermal CFL Higgs phase" 
is dominated by the thermal potential of the instanton size, in which the size is roughly 
equal to the horizon size. The lower half denoted as "CFL Higgs phase" is dominated by 
the one generated by the backreaction due to the baryon density, in which the instanton 
size is much larger than the horizon size. 



new CS-type potential ( p.24| ) adds up on it. This addition does not change the result 
that the instanton size is nonzero, so we still have the Higgs CFL phase. 

If we naively adds up the two potentials (the thermal potential (|5.23| ) given in 
Ref. 1^ and our CS-type potential (|5.24|) ), we find that there are two CFL phases: 
for small values of the baryon density, we have the minimum at a finite value of the 
instanton size, while for large values of the density the minimum sits at the infinite 
size modulus. An expected phase diagram is given in Fig. |^. However, since the 
potential computed in this paper is valid only around small size of the instanton, 
this conclusion is a qualitative one and deserves further study for full backreaction 
of the geometry. 

We conclude by discussing some interesting future directions in Sec. |. 



2. Instanton on the Flavor Branes 

We start with constructing a solution for the equations of motion on the Nf flavor 
D7-branes, which has nonzero instanton number. This solution corresponds to the 
dissolved D3-branes in the flavor D7-branes. In this section, the probe approximation 
for the D7-branes is adopted, while the important backreaction will be treated in 
Sec. ^, and its effect on the solution which we will find in this section will be studied 
in Sec. § where the dynamical dissolution of the instantons (D3-branes) due to the 
finite baryon density is shown. 
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2.1 Review of the D3D7 System at Finite Baryon Density 

Let us begin for simplicity, by considering the case with zero temperature, which 
corresponds to AdS^ x 5*^ background in type IIB Supergravity. We shall embed in it 
a stack of Np space-time filling D7-branes, with a non-trivial world volume baryonic 
f/fc(l) gauge field turned on. As it turns out, the exact shape of D7s and profile of 
the gauge field can be analytically solved in such regime [|14|, which we shall review 
in some detail next. 

The AdSr, x metric, as generated by the backreaction of Nc D3-branes, is 
given in Poincare coordinates: 

ds^ = -^n^udx'^dx'' + ^ {drl + rldsf) , — = Ang^N, = X , (2.1) 

g^d = -^dx'^ A dx^ A dx"^ A dx^ , (2.2) 

gsF^ = (1 + *io)d{gsCA) = iR\dn, + ^lodQ^) • (2.3) 

Here we have listed out the RR 4-form field C4 and the self-dual 5-form field strength 
F5, whereas the string coupling g^ = e'^° remains constant. The indices /i u runs over 
0, 1, 2, 3, 7]fj^^ denotes four dimensional Minkowski metric, and ds^ is the metric for a 
unit five-sphere. For our later purpose, let us also reparametrize the fiat six internal 
dimensional metric as: 

drl + rldsl = dr'^ + r^c/sg + dy"^ + dz^ , (2.4) 

with r| = r"^ + + . Here ds\ {ds^ is the metric on the unit 5*^ {S^)- In such 
coordinates, there exists U{\) C 5*0(6) isometry group which rotates {y,z). 

Introducing Nf (<^ N^) probe D7-branes into (|2.1| ), their Dirac-Born-Infeld 



(DBI) action is given by [|T3] 

SEli = -^D7 j A e-" tT,/-det{Gab + 27ca'Fab). (2.5) 

Here a = 0, ... 7 are the eight dimensional D7-brane world volume coordinates. 
Gab is the pullback metric and Fab is the worldvolume gauge field, which for now, we 
shall only turn on the diagonal baryonic Ub{l) component. TdtC"*" = l/((27r)^a'^(yfs) 
is the tension of the D7-brane. The trace is taken over the symmetrized gauge indices. 
Note that the symmetrized trace is valid only up to the fourth order in a' [0, |l^ . 
However, it is known that, at this order, the non-abelian DBI equations are solved 
at least for the instanton configurations. 

We choose the gauge for the D7-brane worldvolume coordinates as 

(e°,...a = (t,...,x^), {e,...,e)^ir,s'), (2.6) 

so that the D7-branes are spacetime filling and spanning in the four fiat internal 
directions given by r and in (|2.4|) . The D7-branes therefore have asymptotic 
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worldvolume geometry of AdS^ x S^. The precise D7-brane embedding are specified 
by the transverse coordinates {y, z), which become D7-brane scalar fields. To preserve 
the isometry of S^, we have {y{r), z{r)); the U{1) isometry further sets z{r) = 0. The 
induced D7-brane world volume metric is therefore: 

GabdCde = ^''^ ^j^}''^'^ iv.udx'^dxn + , , f, . ((1 + {y'{r)f)dr' + r'dsl) , 

_|_ y[r)'^) 

(2.7) 

where ' denotes Turning on only the temporal component of the ^4(1) gauge field 
Atijr), which we again take to be dependent purely on r, the resultant D7-brane DBI 
action (density)^ is: 

-^^31/^4 = j drL = -U j dr rV(l + (l/'W')-(2v^«'^'^W)^ (2-8) 

where N = NfTj:,i'Vo\{S^)g~^ = NfTQj{2Ti'^)g~^, and factor Nf arises from the trace. 
As noted in Ref. [jl4|, the action (|2.8| ) does not contain explicit dependences on 



y{r) and At{r), their equations of motion yield following constants of motion: 

^ = -ATr^ ^ y' = -c , (2.9) 

Sy' VI + (y'y - i2na'A[y 

\fr^ ^ * — = d . (2.10) 



6{2iTa'A[) + (y')2 _ (27Ta'A[y 

A useful relation can also be readily deduced 

2WA;(r) = 4'(r) . (2.11) 

Using this and rearranging ( p.9|) and ( p.lOj ), we obtain 

2WA;(r) = y\r) = , (2.12) 

where we have defined: 

4 = ^^. (2.13) 



We can readily integrate (|2.12|) to obtain the profiles for y{r) and 2TTa'At{r): 



c 



2 + Vs\ „ ,,,, d ^f,,2 + VS^ 



y(r) = \<^(r), — , 27Ta'At(r) = ^ , \ <^(r,, 

' 2 SVWrg A j ' ' 2 3i/Wr2 " 4 

(2.14) 



^We divide out the infinity volume of four Minkowski spacetime V4. 
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where ¥{ip, k) is the incomplete eUiptic integral of the first kind. In the above 
computations, we have taken d > c, the resultant solutions ( p. 14 ) should be re- 



garded as the zero temperature analog of the aforementioned black hole embedding 



1^, |TJ]. In such case, the D7-branes extend all the way to the "horizon" located at 
_|_ y(7-)2 — — Q (^z{r) has been set to zero) and we have used this fact to fix 
the integration constant. The profile of y{r) in ( |2.14| ) in {r,y{r)) plane displays a 
sharp peak towards y{0) = around r = (or four dimensional cone when sweeping 
out the S^), and flattens out to approach 27ra'm as r ^ oo, where m is the bare 
quark mass. This is in contrast with the Minkowski embedding where D7-branes lie 
at finite distance from the horizon, or ^/r^ + y{r)'^ > th- In the presence of finite 
baryon density, it was shown in Ref. ||13| that only black hole embedding is stable 



and physical, we shall discuss them in more details in section ^ 

Finally one can relate the asymptotic values of y{r) and 2TTa'At{r) with the quark 
mass m and the chemical potential n as y{oo) 2'Ka'm and 2TTa'At{oo) — > 27ia'fi 



[0, and obtain the following relations flj 

c = -f^^{2na'f{^x^ -m^)m, (2.16) 
d = 77V(2W)^(/i^ -m2)/i. (2.17) 

Here the constant 7 = (^ r(i/3^(7/6) ) ~ 0.363. This completes our review on the 
zero temperature D7-brane embedding in the presence of baryonic Uh{l) gauge field. 
To realize the color-flavor locking phase, we shall next consider turning on a SU{Nf) 
instanton configuration within the internal four cycle as a perturbation. 

2.2 Instanton Solution 

We are ready to consider the non-Abelian part of the U{Nf) gauge group on the 
flavor D7-branes, including the instantons. In the equations of motion, the overall 
Ub{l) discussed earlier is coupled to the SU{Nf) subsector where we like to put the 
instantons representing the D3-branes. 

As we shall see later, for large 'tHooft coupling A, the non-Abelian part can be 
regarded as a fluctuation around the fixed Ub{l) background (|2.14|) . We substitute the 
Ub{l) solution ( |2.12| ) of Ref. []14| into the action and consider only the SU{Nf) non- 



Abelian part of the action ( |2.5| ). We are interested in solutions having the instanton 
charges in the subspace (^*^, ■ ■ ■ , ^''), so we just turn on SU (Nf) Ai{^) (i = 4, ■ ■ ■ ,7) 
among the gauge fields, and let them be dependent on only the coordinates i^*. 

In the action, the effective four cycle metric felt by these non-Abelian components 
is computed as follows. We note that the background At{r) can be regarded as an 
additional transverse scalar field in the D7-brane DBI action, as we are interested in 
only the space spanned by (^^, ■ ■ ■ , ^^). Indeed, the effective metric for the directions 
(C^) ■ ■ ■ ; ^'') can be written formally as 

G\f = Qr, + 9vy^^y^,y + g''d,Atd,At{2'Ka'f , = 4, 5, 6, 7. (2.18) 
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Since y{r) and 27Ta'At{r) are functions of r = y X]I=4(^*)^; rewrite above as 

4'^ = f (^^. + ^(^"-(2-«'^'.n)- (2.19) 

So the determinant in the DBI action, including the non-Abehan field strength Fij 
in the SU{Nf), is written as 

- det(G',fe + 27Ta'Fab) = det (d[f + 27Ta' F,,'-^^^'^ (2.20) 

where the unwarped effective four cycle metric clf is given by 

Gff^5,, + ^-^iy"-i2na'A[r). (2.21) 
Thus the total DBI action including the non-Abelian field strengths Fij is 



SEm = -Tm I d'x j d'i e"* tr^det {cf^ + 2na'Fj-^^^^ . (2.22) 

In this expression, note that the prefactor of Fij is suppressed by A~^/^. In fact, 
2'Ka'/R^ = 271 with the relation R'^ = AngsN^a'^. We can therefore regard the 
instanton as a fiuctuation around the fixed Ub{l) background sourced by At, for a 
large A. 

In addition to this DBI action, now we also have a Chern-Simons (CS) term by 
coupling with background RR 4-form C4 given in ( ^.2[ ) with the non-Abelian field 
strength Fij: 

SEl = /^Dr / d'x I d'i 1 (^^^^^) ' tr [e^^^'F^^Fu] , (2.23) 

with /iD7 = We will show that self-dual configurations of the non-Abehan gauge 
fields with respect to the metric satisfies a particular property: the Fj,-dependent 



part of the DBI action (|2.22| ) is completely canceled by the Chern-Simons term (|2.23| ). 



This interesting property of the instantons on the D7-branes was explicitly shown for 



a special case in Ref. ||33| which treated the case of the fiat D7-branes (c = d = 0). 



We use the following formula in generic curved space [36 



^/det{g + F) = ^/d^ + i Vd^ I Fij *4 F'^ \ (2.24) 
for the self-dual configuration 

F^, = *4Fi, . (2.25) 
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Here the Hodge dual operation *4 is with respect to the effective four cycle metric 
and defined by a covariant totally antisymmetric tensor r/*-?'^', 

*F'^ = -Tf^^^Fki , 7]'^^^ = J: e^^'^' , e^^^^ = 1 . (2.26) 

2 ^/det g 



This formula was shown in Ref. for Abelian field strength, and now if we assume 
that the non-Abelian DBI action is written with the symmetric trace prescription, 
this equality also holds for the present non-Abelian case. Once we apply this formula 



to our DBI action ( 2.22 ), for the self-dual instanton configuration with respect to 



the metric G^f , we obtain 



i?2 



ij^ kl 



(2.27) 



Note that we rewrite the Hodge dual by the constant tensor e^^^K Using the relation 
Td7 = /W? and gs = e* is fixed, it is obvious that the FF dependent terms in the DBI 
is canceled by the CS actions ( p.23| ). 

It is interesting that this cancellation occurs not only for the flat D7-branes 
with no electric flux on it but also our present case, albeit our D7-brane configu- 
ration breaks the supersymmetries completely. The state with the D3 branes and 
the D7 branes connected by the fundamental strings in fiat space is supersymmetric. 
However, in our case, the spike does not extend to infinity, supersymmetry is thus 



broken. In Ref. ||33[ , it was argued that this cancellation is due to the BPS property 
of the D3D7 system. Here we could show the same cancellation even with the non- 
super symmetric electric flux, and there is therefore no potential on the instanton 
moduli space. 

However, in the remaining part of this paper, we will see that in fact a back- 
reaction of this electric flux on the D7-branes will lift the cancellation slightly, and 
induces a potential term for the instanton moduli space. It is an essential point 
which we like to focus on in this paper. 

2.3 Conformal Metric and Explicit Instanton Configuration 

Our self-dual configuration of the non-Abelian field strength is with respect to the 
curved "effective" metric G[f . On the other hand, the simpler case of Ref. has 
a fiat metric 6ij instead. In the following, we show that a coordinate transformation 
can turn the effective unwarped four cycle metric G^j^ ( p.21| ) into a conformally fiat 
metric, so that in the new coordinate the standard BPST instanton configuration 
suffices. In any conformally flat space, the self-dual equation on it is simply the same 
as the self-dual equation on the flat space. 



It is easy to see that the metric G[f ( |2.21|) can be written in the polar coordinate 



as 



ds^ = (1 + (y'2 - (27ra'A;)')) dr^ + r^dsl . (2.28) 
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Substituting the explicit expressions for y'{r) and 2TTa'A[{r) (|2.12|) , we can deduce 
that 



l + y'\r)-{2na'A[ir)f = 
To show (|2.28|) is conformally flat, let us consider 



J.6 _|_ 



^6 _|_ ^6 



dr^ + r^dsl = S{r)^ {dr^ + fV,) 



(2.29) 



(2.30) 



and solve for r and S{r). First the consistency in directions demands that r 
S{r)r, the relevant differential equation in the r and f directions then gives: 



dr 



+ 



-dr . 



(2.31) 



Integrating both sides, we can obtain the desired change of variable, 

1/3 



1 + v^l + r{]/r6 



(2.32) 



The integration constant is fixed so that r ~ r for large r. We can also invert the 
relation ( |2.32| ) to obtain 



S{r) 



1 - 



4r6 



1/3 



(2.33) 



In this new coordinate f, the self-dual configuration is just the familiar BPST in- 
stanton. When bringing that to the original coordinate r, we obtain a solution to 
the self-dual equation in the space with the metric G^f . In Sec. ||, we shall use this 
explicit coordinate transformation to evaluate the potential for the instanton size 
moduli. 



3. Linearized Supergravity Backreaction 

In this section, we shall compute a linearized perturbation to the supergravity 
background ( |2.1|) , ( p.2|) , (|2.3| ), due to the electric field At on the D7-branes. Let us 
first recall that the electric flux, which is responsible for the Ub{l) baryon charge, 
can be regarded as fundamental strings dissolved in the D7-branes. This is because 
in the DBI action the electric field is combined with the (pull-back of) NSNS 2-form 
field i?2 in a gauge-invariant fashon, 27ca'Fab + Bab- Such electrified D7-branes can be 
regarded as a source to the bulk 3-form flux = dB2, acting as small perturbation 
to the background SUGRA solution. Moreover from the consistent equations of 
motion of the SUGRA, this also induces RR 3-form flux F^, which we will proceed to 
extract in two different ways. The induced is important for the dynamics of the 
instantons on the D7-branes as we will see in the next section. So, in this section. 
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we derive the exact amount of this as a backreaction of the electrified D7-brane 
configuration, which is 

Fili-'-^- (3.1) 

First in Sec. we present an intuitive derivation of the F3 by using smeared 



baryon vertices. In Sec. |3.2| , we compute the backreaction to the geometry due to 



the electrified D7-branes. The result of Sec. 3.2 coincides with that of Sec. 3.1 



3.1 Smeared Baryon Vertices 

The electric fields on the D7-branes are interpreted as fundamental strings connecting 
the D7-branes and the D3-branes, therefore they are quarks. The number density of 
them is given by ( p.lO| ), quark density = 2iTa'd. This means that the baryon number 
density is 27ia'd/Nc. 

The D7-brane spike terminates at the origin tq = 0. If we take the flux con- 
servation at the tip of the spike seriously, we need to assume the presence of the 
baryon vertices surrounding the origin. As is well known, D5-branes wrapping the 
S^, which are called baryon vertices, can give a charge at which the fundamental 
strings can end [Q. In this subsection, we compute a back reaction of these baryon 
vertices smeared on the plane x^-x'^-x^ at rg = 0. Our result is (|3.1|) .^ 

The relevant terms from the type IIB supergravity and D5-brane DBI actions 

are 



4k2 



d'"x ^^o\F7\' + 1^, / a, (3.2) 



with = 2(27r)^a'^ and fi^ = (27r)"^a'"^. We have also used F7 = dC^ = *ioF3. It 
is enough to consider the explicit component Cg = C^0^g^0,_ig^g^dx^ A dOi A d62 A dO^ A 
dO^ A dO^, then (B.2p becomes 



d'^xdrQd6id02d0^d0iid0r^ 



-1 r| 1 



4Kfo sin^ ^1 sin^ 6^ sin^ ^3 sin ^4 ^'"''^'^oe.e.eae^e. 



"10 

2Wd 



(3.3) 



where the position of the baryon vertices is specified as rg = e with e 0. This can 
be solved as 

^reCSe^eae^e^ = — -3- sin^ ^1 sin^ 62 sin^ ^3 sin 6^ , (3.4) 

where we have also used the explicit expressions for k\q and /is. Taking a Hodge 
dual in the background AdS:^ x S"^, we immediately obtain (|3.1| ). 



related issue on backreaction of baryon vertices was discussed in Ref. [Eol 
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The above analysis leads to an important consequence which resolves a problem 
in introducing baryons in Dp/Dg systems. The phase structure of fundamental mat- 
ter at finite baryon density has been studied by introducing electric flux on probe 



Dg-branes in Dp-brane background [1^, [T^, [S^] • The baryon number there was con- 
sidered to be carried by free quarks in the sense that the quark density, or the electric 
flux, can take any value as long as the total number of strings takes an integer. In 
other words d is quantized in units of 1.^ It is natural to ask what if one considers 
baryons instead of the quarks in the system. As was pointed out in Refs. |T^, ^ 
and studied in detail in Ref . , it turned out that there is no stable baryon vertex 
solution outside the horizon, in the deconfinment phases. 

A resolution of this problem of the missing baryon vertex is that the baryon 
vertices undergo a brane/flux transition and leaves only RR flux outside the horizon. 
The DBI part of the D5-brane baryon vertex disappears since its volume element van- 
ishes (the time direction of the geometry shrinks), while the CS term of the D5-brane 
action remains to source the bulk RR 3-form flux F3. We can see this "remnant" of 
the baryon vertices anyway, by solving consistently the SUGRA equations of motion 
for the NSNS -B-field, as in the following Sec. |3.2| . In this section we work with 
T = 0, but the role of the horizon is played by the origin = 0. Similar transition 
can be found in a simpler example. Consider AdS^ x with N units of the RR flux 
and put an additional probe D3 brane parallel to the boundary in this spacetime at 
certain r^. This is a supersymmetric configuration (the Coulomb phase) and rg is a 
modulus. When the brane goes to rg = 0, the DBI part of the D3 brane becomes 
zero. The correct picture for this case is given by AdS^ x 5*^ with + 1 units of the 
flux. Therefore, the probe D3 brane is replaced by a unit of flux. This argument can 
be applied to the finite temperature case. 

It is interesting that the quantization condition of this -F3 in ( p.l|) shows that 
the quark number density d is quantized not in units of 1 but in units of Nc 
This would suggest that the quarks in the Dp/Dq system are always thought to be 
components of baryons. We will see in Sec. |^ that the analysis of Sec. still applies 
to a finite temperature system despite the fact that the end points of the strings are 
hidden inside the black hole horizon. 



3.2 Backreaction from the D7-brane Electric Flux 

Instead of assuming the presence of the D5-brane baryon vertices, here we provide 
an alterative derivation for (|3.1|) by solving the backreaction due to the electric flux on 

^In this paper this d is the density, but one can imagine locahzed quarks/baryons instead, for 
the discussion here. 

^"The flux is smeared along the space directions parallel to the boundary: xi,X2,X3. With these 
directions being non-compact, we do not need to quantize the flux from computational point of 
view. However, our motivation of this quantization comes from the fact that the flux is sourced by 
the D5 branes. 
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the D7-branes. In this subsection, we demonstrate this by looking at the equation of 
motion for the NSNS 2-form field B2. For the validity of our approximation adopted 
in this section, see Sec. 0. 

3.2.1 Sourcing the Bulk NS-NS B-Field 

First, let us examine how the electric flux on the D7-branes can act as a source for 
the bulk NSNS 2-form field B^. The DBI action includes the NSNS B-field as 



qD7 
'-'DBI 



Tor / rf^e e-* tiJ-detigab + 27Ta'Fab + B, 



'ab, 



(3.5) 



where Bab is the induced NSNS 5-field carried by the fundamental strings in the 
D7-branes. We are treating here only the overall U}j{l) C U{Nf) sector (|2.12| ) only, 
so tracing over already gave the factor Nf. Since we are interested in a linear per- 
turbation by this source, we expand this action around i? = to the linear order in 
the 5-field: 



cD7 
•^DBI 



0{B) 



— / d^xdr B\ 



Or 



5L 



(3.6) 



B=0 



where L is the Lagrangian density as defined in ( p.8|) . This expression follows since 
the S-field appears in the DBI action ( |3.5| ) only as a gauge invariant combination of 



2TTa'F + B. More explicitly they appear as [211 a' A ^ — -Bor)^ in the action, hence we 
have the additional negative sign. The shape of the D7-branes is specified only by 
y{r), so the induced 5-field is just = BQyy'{r) + i?or- Together with ( ^.1(J| ), we 
have 



^DBI 



0{B) 



-d j d^xdr (Boyu' + Bor) ■ 



(3.7) 



This is the source term for the bulk NSNS B-field. 

For our later purpose, we can express (r, y, z) in terms of angular coordinates for 
the ^5 {^i,...,^5}: 



Tq cos 61 , y = Tq sin 6*1 cos 6*2 , r = rg sin 61 sin 6*2 



(3.^ 



the remaining 63,64,, 65 parametrize in ( |2.7| ). For our D7-brane embedding spec- 
ified by y{r), z = 0, this translates into setting 61 = 7r/2. We can further invert the 
relation rg = + y(r)^, and express 62 as a function of rg via: 

r(r6) 



62{fe) = arctan 



y{r{re)) ' 



(3.9) 



So in terms of rg and these angular coordinates, the source coupling for the NSNS 
5-field is 



S 



D7 
DBI 



OiB) 



-d I d X I drQ diZs ^ „ . 3^ ^org + Boe^jr- 

27r^sm'^6'2 V '^'"e. 



.(3.10) 
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Note that to incorporate the D7-brane DBI action into the full 10 dimensional su- 
pergravity analysis, we have inserted delta-functions which restrict to the specific 
embedding we are considering. In particular we have changed the integral to the 
whole angular coordinates, so we divide it by the volume V3 = 27r^ of the unit 
3- sphere. 

We will now concentrate on the region r ~ to simplify our situation. The spike 
has a rigid cone shape around the origin r = 0. Around the tip of the cone r ~ 0, 
diverges, so around there we have 



With this, we can approximate the source term ( p.lOp as 



cD7 
^DBI 



0{B) 



d^x / rfre cifis - 7^/2)5(^2 - 0. 



(0)^ 



5, 



Ore 



sin^ef ■ 



(3.11) 



(3.12) 



3.2.2 Extracting the RR 3-Form Flux 

We now would like to extract the linearized perturbation F3, which will be crucial for 
generating the potential on the instanton moduli space. For this, at the linear order 
perturbation, it is sufficient to consider the equations of motion for the NSNS 5-field, 
with this limiting source term ( |3.12| ) included, as other equations are affected only 
at higher orders (we will check this later, see eq. ( |A.1| ) and discussions thereafter). 
It will be important that near the r ~ region, this source is only for -Bore- 
The relevant part of the type IIB supergravity action is |3^ 



^'^lO 



Here 



so, in the third term in ( |3.13| ), the term linear in B2 is 

^ ' ' -C2AH3 + B2A F3) A *F. 



10 



5 • 



(3.13) 



(3.14) 



(3.15) 



*F. 



this is equal to the second 



For a self-dual background 5-form flux F5 
term of (|3.13|) . 

Substituting the AdS^ x background metric and the RR 5-form flux (|2.3| ) and 
writing out in explicit components, we obtain 



1 



2(27r)V4(?2 
1 



d^^xdr^dVL^ rl 



■2 n ^-"-Oreei 

sm (79 



+ 



,/4 



d^xdredQ^ BQr^F^^2'^T:N^{af . 



(3.16) 
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Here we used the explicit 5-form flux on the S^, = 2'^Ts:N(.a''^dVL^ (|2.3|) . Together 
with the source action ( p. 12 ), the total equation of motion for the NSNS S-field is 



= ^7' M 2 [^^1 (sin^^isin^^s^^oreej +5^. (sin^ sin^ ^^a^^oree,)] 

-Fif^ sin^ 01 sin^ 02 2^T^N^a'^ 

d 

'272 



"(27r)V4^23' 

-5(^i-7r/2)5(^^2-^f )A- (3-17) 



We can recognize this as a 1+2-dimensional electromagnetism on a compact space 
spanned by 9i and ^2- The first term is a total divergence, so the remaining terms 
should vanish when we perform an integration over the 2-dimensional space. This 
condition results in 

TTTTT^^S FdO^ sin^^i f (162 sin=^ ^2 (2ViV,«'2) = . (3.18) 

Performing the integration and re-arranging, we obtain the constant RR 3-form flux 
F^^23 given in ( |3.1|) . This is the leading order effect of the backreaction of the 
D7-brane electric flux. Note that the supergravity equation of motion for the F3 flux 
is trivially satisfied with this constant configuration. 

Interestingly, this result ( p.l] ) is the same one obtained previously by solving the 
F3 equation of motion with the smeared baryon vertices in Sec. |3.1| . Here we have 
not assumed the presence of the baryon vertices, but the supergravity equation of 
motion "knows" the presence for its consistency. 

4. Dissolution of the Instanton and Color-Flavor Locking 

In this section, we show that the backreacted supergravity flux ( p.l|) generates a 



nontrivial potential for the instanton moduli space (Sec. [4.1|) . It provides a dynam- 
ical mechanism for fattening the size of the instanton on the flavor D7-branes. We 
compute the potential explicitly in Sec. 



This is a dynamical color-flavor locking in the holographic QCD, because the 
size of the instanton on the D7-branes is the vev of the squark of the supersymmetric 
QCD. Since the squarks are in the bi-fundamental representation of the color and 
the flavor symmetries, their condensation gives a color-flavor locking. The squark 
condensation means that the theory favors Higgs phase when the baryon chemical 
potential /i is larger than the quark/squark mass m. We study the patterns of the 
symmetry breaking in Sec. [4.3| . 



4.1 Additional Chern-Simons Term 

Using the backreacted supergravity solution ( p.l|) , we obtain an additional Chern- 
Simons term induced on the D7-branes. The general formula for the Chern-Simons 
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couplings on the D7-branes is 

/i7 tr / exp(27ra'F + B2) A^Cg. (4.1) 

Here the D7-brane RR charge is fi-r = ^^^^^^,4 , and the field strength 271 a' F now 
also contains non-Abelian instanton piece. Formally expanding ([4.1|) out and per- 
forming integration by parts, for non-zero F]^^23 obtain the explicit expression in 
components 

^cs = ^^^^ j d'xF^, j d^i tr [A,F,,F,,e^=^^] + • • • (4.2) 



Here • • • means terms necessary to form a gauge-invariant CS 5-form 

tr 



AFF - -A^F + —A^ 
2 10 



(4.3) 



where the wedge product A is omitted. The second and the third terms in the CS 
action are irrelevant for our subsequent discussions. We substitute the constant F^ 



(3) 
123 

from the linearized supergravity backreaction (3J.) to extract the relevant term from 



a'd 

16^ 



d^x j d\ tr [AoF,,F,,e^^'''] . 



(4.4) 



This is the leading correction term due to the baryon density d. Note that the factor 
1/Nc in front of above shows that this is indeed a correction to the original D7-brane 
action. 

This additional CS term has an important physical meaning. The essence here 
is quite similar to the generation of the baryon charge in the Sakai-Sugimoto model 
[0, while the use of the instantons here is quite different from there (c. /. footnote 
6.). As discussed in Sec. |l|, we are studying the process of moving one D3-brane 
from the origin onto the worldvolume of the D7-branes. Once the single D3-brane 
goes outside the D5-brane wrapping the 5*^, the RR 5-form fiux penetrating the 5*^ 
worldvolume of the D5-branes reduces by one unit. This 5-form was responsible for 
the CS term on the D5-branes to produce the electric charges on the D5-branes, 
which are the end points of the fundamental strings. So, by this moving process, 
the total number of the fundamental strings decreases by a fraction of l/N^. Then, 
where does the baryon charge go off to? The answer is the new CS term ([4.2|) . Once 
the D3-brane gets on the D7-branes, it induces an instanton charge. The instanton 
number for the single instanton is 

tr / d^i F,,Fue'^^' = 'i2'n\ (4.5) 
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Figure 4: A schematic picture of tlie cliarge conservation process. In the left figure, 
D3-branes are sitting inside the baryon vertex (D5-brane wrapping S^). From the D3- 
branes, Nc units of RR 5-form flux emanates. The CS term on the D5-brane creates an 
electric charge on the worldvolume of the D5-brane Q , and this generates electric field on 
the spiky D7-brane which touches the D5-brane. In the right figure, we move one D3-brane 
toward outside of the baryon vertex. This D3-brane becomes an instanton (shaded region 
on the D7-brane spike). The instanton is electrically charged, so the total electric flux 
going to the asymptotic infinity of the D7-brane worldvolume is conserved. 



d^xd""^ tiAt 6\0 ■ (4.6) 



For a small size instanton, the CS term ( |4.2| ) effectively becomes proportional to 

27ra'd 

This means that the instanton (which is the D3-brane dissolved into the D7-brane) 
carries the electric charge {2na'd)/Nc. Compared to this amount of the charge, 
the original solution ( p.l2| ) shows that the spike has an electric charge 27Ta'd. We 
can therefore conclude that, pulling out one from the Nc D3-branes decreases the 
baryon charge by its fraction 1/Nc. This decrease is indeed offset by the instanton 
sourcing the electric field, so that the asymptotic expression for the electric field 
doesn't change. See Fig. ^ for a schematic explanation of this charge conservation. 

4.2 Potential on the Instanton Moduli Space 

Finally we have collected all the pieces for computing the induced potential on the 
instanton moduli space, and we will show that it drives the instanton(s) into disso- 
lution. Here we consider a simpler special case, where the BPST instanton profile is 
centered at the origin: 

trF.,F,,e-- = j^^^ . (4.7) 

Note that the BPST instanton solution is obtained in the conformally flat r coordi- 
nate, not in the r coordinate. Upon substituting into the additional CS term 
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we can extract the potential for the size p, Vb{p), via the relation S = — J d'^xV{p) 



as 



12p'(Wd) f3 



vM^- 'Kr^' I *AM-^. (4. 



Again, note that the argument of the electric potential At is r which is related to r 
by (|2.32|) . Integrating by parts (for the f coordinate) and use (|2.12| ), we obtain 



Vb{p) = ^ dr A;(r 



= _2Wd d p^(3r-^(r)+p^) 

A^AT, X (r~2(r)+p2)3 • ^ 

This (^4.9|) is indeed a monotonically decreasing function of p, when viewing together 
with the coordinate redefinition (|2.32| ). This can be easily understood if we notice 
following three facts: (i) A'^{r) is a monotonically decreasing function of r. (ii) The 
last factor in the integrand of ( [4 .91) is the instanton density function which is peaked 
at r = and monotonically decreasing in f, while the width of the function is given 
by p and the function has a normalized integral (which is the instanton number), 
(iii) The map ( p^.32[ ) between r and f is a one-to-one and monotonic function. The 
instanton modulus potential V{p) ( [4.9|) shows that the system dynamically favors 
the Higgs phase, p 7^ 0. This is the color-flavor locking in the supersymmetric QCD. 

The Vb{p) from the CS term computation ( [1.9|) does not appear to be analytically 
integrable, however to get a qualitative understanding we can consider the following 
asymptotic values: 

Vb(p = 0) = 0, Vb{p = ^) = ^/ Kdr = (4.10) 

Their difference, 



Nc Jo 



Vb{p = 0)-Vb{p = oo) = ^-^ (4.11) 

is consistent with the interpretation that we pull out one quark per each baryon to the 
infinity in the background chemical potential p. However, note that the constant -F^23 
is obtained only in the vicinity of r = 0. So, our calculation is strictly valid only for 
small p, and the potential height ( [4.11| ) derived at large p should not be reliable. We 
however expect that the qualitative physical result ( [4.11[ ) is not modified significantly 

(3) 

when we taken into account of full at large radius. For completeness, here we 
include the plot of the one instanton size modulus potential, normalized by the 
asymptotic value Vb{oo): In the analysis above, we assumed that the center of the 
BPST instanton is at the origin. However, we can consider a generic situation where 
the center of the instanton is not at the origin r = 0. Suppose that the center is 
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at some distance X from the origin. It is clear that the similar expression to ([4.9| ), 
which is again only valid at small r, would give a qualitative result that the potential 
Vb(/5, X) goes to its minimum at (p, X) = (cx), finite) or (p, X) = (finite, oo). The 
latter is in particular an extreme point in the moduli space in the Coulomb phase. 
There p can vanish and in that case the color-fiavor locking does not occur. However, 
for X 7^ 0, the original gauge group U{Nc) is explicitly broken, as it is in the Coulomb 
phase. 

Instead of substituting the BPST instanton, we can substitute multi-instanton 
solutions. Suppose we treat Nc instantons. We need to require Nc -C Nc, because in 
the perturbation of the backreaction the original background of AdS^ x should not 
be drastically modified. Generically, the instantons are separated from each other. 
Our analysis for the single BPST instanton holds also for the multi-instanton case. 
Then it is shown that all the size moduli of the instantons are destabilized. 



4.3 What is Locked? 

Naively speaking, the condensation of the squark field which comes from a fiuctuation 
of a string connecting the D3-branes and the D7-branes gives the instanton size. This 
string transforms in the fundamental representation of the gauge group U{Nc) and 
in the anti-fundamental representation of the fiavor group^^ U{Nf). Here U{Nc) is 
a part of the total gauge group U{Nc) i.e. Nc < N^. The partial color symmetry is 
for our convenience to consider instantons only to dissolve into the worldvolume 

"'^'^Since electrically charged instantons should repel each other, generic configurations should be 
with the separated instantons. This phenomenon is common with the interaction among baryons 
in the Sakai-Sugimoto model the repulsive core of nucleons is mainly due to this electric 
repulsion. 

-'^^Precisely speaking, the fiavor group of the supersymmetric QCD should be SU{Nf), since the 
overall U{1) transformation of the fiavor symmetry can be identified as a global part of a U{1) 
subgroup of the local gauge group. In the following, we adopt U{Nf) rather than SU{Nf), but the 
argument goes similarly. 
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of the D7-branes. By the squark condensation, apparently the color and the flavor 
symmetries U{Nc) x U{Nf) are broken. We shall see the symmetry breaking pattern. 
The theory on the instantonic D3-branes is constructed from their ADHM data 
The ADHM data consists of four U {N^) adjoint scalar fields which are combined 
into two Nc X complex scalar fields Bi and B2, and the squark fields P and J 
which are complex Nf x Nc matrix scalar fields. The squark fields transform as 

P ^ UPUq\ J ^ UJUq\ (4.12) 

where U G U{Nf) and Uq G U{Nc)- So these fields are in the bi-fundamental 
representation. 

Let us consider 't Hooft instantons. Then Bi and B2 are diagonal matrices, 
and the ADHM equations are nothing but the BPS equations for the theory on the 
D3-branes, are 

lP = fj, IJ = 0. (4.13) 

This of course allows a trivial solution P = J = 0, which corresponds to the zero-size 
instanton. 

First we consider the two-flavor case Nf = 2. The flavor symmetry is U{2) which 
is a vector part of the chiral symmetry (the chiral symmetry is explicitly broken by 
the quark/squark mass in our case). The 'tHooft instanton whose centers are located 
at the origin is represented by a solution 

Here, all p^'s are real parameters. These correspond to size of each instanton. With 
this at hand, we can compute unbroken symmetry. 

For simplicity, we consider the case of a single instanton, A^c = 1. 



.0; VI. 

Here p is a nonzero constant (for which we computed the potential). In this case, 
the transformation which leaves P intact is 

U={ ^ ,„J, Uo = e'''\ a. gM. (4.16) 



^ e*"2^ 

On the other hand, the symmetry which leaves J intact is 



U=[ ^ ), ^^o = e*"^ a. gM. (4.17) 



e 
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Therefore, we need to require ai = 0^2, and we conclude, 



f/(l)color X f/(2) 



flavor 



CFL 



(4.18) 



The global part of t/(l)coior is locked with the diagonal part of t/(2)flavor, which is 
the baryonic symmetry U{1)b, and the local part of f/(l)coior is broken. 

This kind of the color-flavor locking can be found for general Nf. In the case of 
Nc = 1, the squark condensations are given by 



Jt 



P 



/1\ 







J 



1 





(4.19) 



A similar analysis shows 

f/(l)color X f/(Ar^)flavor ^ f/(l)cFL X U{Nf - 2)flavor • 



(4.20) 



^(l)cFL is a global symmetry which locks a part of the flavor symmetry and the 
gauge symmetry. Note that for Nf > 2, this f/(l)cFL can be nothing to do with the 
baryonic symmetry U{1)b, since the action of this f/(l)cFL can be chosen as 



U = diag(e*",e*",e- 



-2ia/Nf 



-2ia/Nf^ 



Uo 



a G 



(4.21) 



For generic N^, we expect that all the size moduli are driven to have nonzero 
values (for example, for Nf = 2 we have ( [4.14|) ). So the symmetry is broken as 



t/(iVj,olor X U{Nf) 



flavor 



f/(l)cFL X U{Nf - 2) 



flavor • 



(4.22) 



The locking is quite similar to the case of A^c = 1- Note that we restrict ourselves 
to the case of the 'tHooft instantons. Since the 'tHooft instantons do not cover all 
the moduli space of the instantons, there remains a possibility that the unbroken 
symmetry, in particular the part concerning the gauge symmetry, may be enhanced. 

In the D-brane analysis, we used the technique for Coulomb branch in AdS / CFT 
correspondence and treat one D3-brane as a probe by separating it from the rest, by 
hand. This procedure for the separation is somewhat artificial, but it is required for 
the geometry not to be drastically modified by a possible back reaction, which is our 
limitation. 



5. Extension to Finite Temperature System 

In this section we explore the system at finite temperature and baryon density. Since 
the boundary geometry of our AdS^ x is i?^'^, the geometry creates a horizon 
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inside AdS at any finite temperature. A/" = 4 SYM theory coupled with M = 2 
hypermultiplets in fundamental representation at finite temperature has been studied 
(See for a review). As mentioned in Sec. there are two brane embeddings in 



the black hole background: the Minkowski embedding and the black hole embedding 
(Fig. 0). The former describes the configuration of D-branes staying outside of the 
horizon everywhere and the latter describes the configuration of D-branes falling 
into the black hole horizon. The equations of motion and the free energy determine 
which configuration is realized at given quark mass and baryon density normalized 
by temperature: m/T and 27ra;'d/T^. 

In the case of zero baryon density and no instanton, the Minkowski embedding 
covers only higher m/T region while the black hole embedding covers only lower 



m/T region ||26|. These two are connected by a first order phase transition at cer- 
tain critical temperature (m/T) crit.- This phase transition is interpreted in the field 
theory as meson melting: the spectrum is discrete and the mesons are stable in the 
Minkowski embedding, while the spectrum is continuous and the mesons are unsta- 
ble in the black hole embedding. In the case of finite baryon density, the Minkowski 
embedding is no longer physically allowed and the black hole embedding covers the 
whole temperature region. When an instanton is excited on the D7-branes at zero 
baryon density, the potential for the instanton size moduli takes its minimum at the 
origin, p = 0, for the Minkowski embeddings and at some finite value, p = p^am > 0, 
for the black hole embeddings This means the system is in Higgs phase above 
the critical temperature (m/T)crit.- Therefore, the system is already in a CFL phase 
with the squark condensation, in the melted meson phase with finite T. 

The purpose of this section is to study the case with both the finite baryon 
density and the instanton configuration. As we saw in the zero temperature case, 
the backreaction of the baryon density excites an additional CS term which induces 
the CFL. We will see how this CS term affects the instanton potential in the finite 
temperature system. 

5.1 The D3D7 System at Finite Baryon Density and Temperature 

The background geometry dual to the finite temperature system is an AdS black 
hole. In Poincare like coordinates, the metric, RR 4-form and the dilaton have the 
following expressions, in the conventions of Ref. |TB[: 



rfs^ = (^-Ldt^ + fdxl^ + ^ {du^ + u^dsl) 

2 4 \ 2 



where 



/=l-4. /=l + 4. (5-2) 



- 24 - 



with Mo the location of the horizon, u and Tq in Sec. 3 are related to each other by 
the coordinate transformation 



u' = ri + ^ri-ul (5.3) 

The regularity of the Euclidean section of this geometry relates the horizon radius uq 
and the Hawking temperature, which is interpreted as a temperature of the boundary 
gauge theory of our concern, as 

^ = ^ • (5-4) 

As in the zero temperature case, we foliate with 5''^ so that 50(3) R-symmetry 
is manifest: 

du^ + u^dsl = du^ + u'^ide^ + sin^ Odsj + cos^ 9d(j)'^) , (5.5) 
where 6, and 6i, 62 in ( p.8| ) are related through the following equations 

sin 6' = sin 6'i sin 6'2 , tan0= ^— — . (5.6) 

sin 9i cos 02 

The probe D7-brane worldvolume is spanned by (t, x*, m, S"^), and is localized in 
direction, which we can use rotational symmetry to set = (corresponding to 
9i = |). In the new coordinates, the embedding is described by x = cos 6* as a 
function of u. The Uf,{l) gauge field A dual to the baryon current on the boundary 
has only non-zero time component 

A = At{u)dt. (5.7) 
With these ansatze, the DBI action of the D7-branes is given by 

du C 



cB7 
^DBI 



= -N J du W 1 - X + - {27Ta'Aty — (5.8) 

where the dot denotes ^ and Af is as defined below ( p.8|) . Since the action does not 
contain A explicitly, the momentum conjugate of A is constant: 



or equivalent ly, 

' 7(1 - x^) Jn6/^(l - X^? + 8(D/Ar)2 
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To derive the equation of motion for x, we Legendre transform the action with respect 
to D so that At can be ehminated from the action: 

C = C ^^^(2WA) 

6{27ra'At) ^ ' 

AT / 
Then the x equation is 



1 - + u^x^\Ju^P{l - x^f + 8(D/Ar)2 . (5.11) 



^ 8(D/Ar)^ 



'1 - + V "^^^(1 - ^')^ 

X I 3(1 - x') + 2mV - 24 V . + I _ (5^^2) 

' ^ ^ VA/^y MV=^(l-x')=^ + 8(D/Ar)2 ' ^ ^ 



As studied in Ref. the boundary condition of probe branes at the horizon is de- 
termined by the regularity of the induced curvature: xk=«o — 0- With this boundary 
condition, the solution of this equation of motion near the horizon m ~ mq is then 
given by 

Therefore, the embedding can be approximated as 

X = Xo, X = 0, (5.14) 
^ ~ 1 smaller than ^(i^tg)^/^ Ar?/^ ^^^^ ^^^^e D. 

We again consider the instanton excitations as a perturbation in this background 
field. At the leading order of D/A^^c, the instanton couples to RR 4-form and the 
induced metric, and the relevant terms are: 

S^^iiPP) = -NfTr,, I d'x I ^Jf- Tr[F A F] , (5.15) 

S^l{FF) = Nf%, I d'x I • Tr[F A F] . (5.16) 

The instanton F A F lives on an effective four dimensional space whose metric is 




^ 2 / 
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^ ' -du^ + - x )dsl 



(5.17) 



Note that D dependence appears only through this metric. The difference between 
the DBI term and the CS term is the factors of / and / in the integrands. Therefore, 
the instanton potential vanishes as long as the temperature is zero, even in the 
presence of finite baryon density as we have seen. 

5.2 CS Term from Backreaction 

The next leading order in Ti/Nc comes from the backreaction of the gauge field to 
RR fields. Similar calculations show that it excites the same constant RR 3-form 
field near the horizon as that of the zero temperature case with d replaced by D, 

41^^. (5.18) 



c 



Therefore it induces the same CS term 

= 1^ I rf^a: I tr [A A F A F] . (5.19) 

The remaining issue is to obtain the explicit instanton configuration tr [F A F] on 
the effective metric Gj-j"*. As in the zero temperature case, we would like to obtain a 
conformally fiat coordinate u satisfying: 

Gg) = Siufidu" + u^ds:,) . (5.20) 

For the approximate embedding near the horizon ( |5.14| ), the conformally flat metric 
( p. 201 ) is obtained by 

J^ltzM^=du='±. (5.21) 

With this conformally flat coordinate, the BPST instanton is given by 

tr[FAF] = ^^l|^rf^e. (5.22) 

As we have seen, there are three terms contributing the instanton potential for 
the finite temperature case: the DBI and the CS term for thermal effect, which we 
denote Vr, and the CS term from the backreaction, which we denoted earlier as Vb- 
From ( ^.15 ) and ( |5.16|) , the thermal potential is 

Vt{p) = ~iSEl,{FF) + S^l{FF))/V, 
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= -4^ y - /(/ - /) ■ 8 (^^+py ■ ^ ^ ^ 

This potential has a minimum at finite p since Vt(0) = Vt{oo) = 0. On the other 
hand, from (|5.19|) , the potential from the backreaction is 

Vb{p) = -SO'(D)/14 

, V '';'!!^;?.(5.24) 

iV,A/ JfJu^p^l - ^1)3 + 8(D/Ar)2 +P > 



u=uo 

The approximation we used to obtain the integrands breaks down for large u. 
However, the following two features may still hold even beyond the approximation: 
Vb{0) = and Vb(oo) = _{?ZI^i^^ xhe first one, Vb(0) = 0, comes from the 
fact that the integration of the instanton term, the last term in ( |5.24| ), is zero for 



p = 0. Therefore, independent of the form of F3 and At, the integration gives zero. 
The second one, ^3(00) = _(Hzi^i2/£^ comes from the physical reason explained in 
Section 4.2 for the zero temperature case. A numerical analysis suggests that Vb 
monotonically decreases from zero to _(?Zl^i^_ 

The shape of the total potential V = Vt + Vb then depends on the ratio between 
them, which is characterized by 



4 



^ I \ (5.25) 



Vt 1271 \Af 

Physically, this suggests that the thermal effect dominates when the temperature is 
high (large Mq), while the backreaction effect dominates when the baryon density is 
high (large D).^^ Note that (D/A/") can be large up to the order of {Nc/Nf) where 
the probe flavor brane description breaks down. Therefore, the ratio (Vb/Vt) may 
become large despite the fact that it is a positive power of gg. Recalling that Vr 
has a local stable minimum and Vb is a run-away type potential, we conclude that 
the potential can have three possible behaviors depending on the ratio between D 
and Uq. When ((D / J\f)'^^^ / Uq) is very small, Vt dominates the potential and it has 
a local stable minimum. The size of the instanton in this case is about the order 
of the horizon scale. Since observables which have less energy than the tempera- 
ture have no meaning at flnite temperature, the flniteness of the instanton size may 
be interpreted as a thermal effect. As {(D / J\f)^^^ / Uq) increases, the local minimum 
becomes a meta-stable state and the instanton size p eventually decays to inflnity. 
As {{'D/^^y/^/u^) increases further, Vb dominates the potential and the local min- 
imum disappears. These features are shown in Figure |[ Thermal quantities such 
as the derivative of the free energy may change discontinuously at a critical value 
of {{D / J\f)^^^ / Uq) . Therefore, this shows a phase transition within the CFL phase. 



13 



Note that both potentials arc of 0{1/Nc) for fixed gs (~ fixed ratio of X/Nc)- 
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V(p) 




Figure 6: The plot of i^^^^-^-^i versus for xo = 2 • '^^^ three Unes, from bottom to top, 

correspond to ^^^3 =0.4, 0.8, and 10 respectively. We can see that for when the baryon 
density is small compared to the temperature, the thermal potential Vt dominates and 
the potential has a local stable minimum. As the baryon density is increased, the relative 
contribution of Vb becomes larger and eventually the potential becomes the run-away type. 



For larger temperature, we have a CFL phase with the finite size instanton, while 
for larger baryon density, we have another CFL phase with the size of the instanton 
being very large. A schematic picture of the phase diagram is shown in Fig. |^. 

As in the case of T = 0, our analysis is valid only for small u region and cannot 
say anything precise about the potential at large u. From the boundary theory 
point of view, the stability at higher temperature may be understood as the thermal 
masses of the scalar fields and the instability at higher density may be understood 
as the tachyonic masses of the scalar fields from the chemical potential. Since the 
supersymmetry is completely broken, the potential might be lifted up by cubic and 
higher terms and the vevs of the squarks may take finite values. Of course, these 
expectations are from the weak coupling analysis of the gauge theory and the strong 
coupling dynamics might change the picture. We do not go into the detail on this 
point in this paper. 



6. Discussions 

In this paper, we made some preliminary steps towards a holographic model of color- 
flavor locking phase, here we end with a list of interesting future directions which 
seem worth exploring. 

The phase diagram (Fig. |^) is obtained by the total potential Vr + Vb for the 
instanton size modulus, but the potential Vb is valid only for a restricted region for 
r, as shown in Sec. |^. So, it is important to compute the backreaction which is valid 
in all region of r, to explore the phase diagram further. 

In particular for T = 0, we have shown that there is an instability along the 
direction of squark VEV in the melted meson phase. This means that the critical 
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chemical potential dividing the meson and melted meson phases may take a different 
value which is smaller than fi = m. Our method of treating D3-branes among 
Nc of them separately cannot reach the true value of the critical chemical potential 
in the full phase diagram, and this deserves a further study. It is possible that there 
may be no vacuum if the potential valid for all r is found and turns out to be a 
run-away type. See |^ for a related discussion for R-charge chemical potential. 

A related issue is a possible distinction between the two CFL Higgs phases. We 
have two CFL phases, one is with finite instanton size p while the other is with p = oo. 
The former is realized mainly by the thermal potential for p, while the latter is by a 
domination of the baryon density. The symmetry breaking patterns look similar to 
each other. However, we expect that, once the repulsion among electrically charged 
instantons is included, the remaining symmetries may differ. In addition, physical 
solitonic spectra in these CFL vacua may be different from each other. It would be 
interesting to study vortex strings in these vacua. 

The vortex strings in the CFL phase in QCD play important roles in various 
physics (see Ref. for a partial hst of related papers), and the D-brane techniques 
for the CFL phase studied in this paper may be helpful in revealing the properties 
of those vortex strings. Since the vortex strings are inside the D3-branes which 
are instantons on the D7-branes, this suggests that "vortices inside instantons" are 
possible. This is intriguing on its own in soliton physics. 

It was described in Ref. that in an idealistic situation the CFL phase of QCD 
may be continuously connected to the hadron phase, giving a continuous deforma- 
tion of the excitation spectrum, named "quark-hadron continuity". In our case, the 
dynamically favored CFL phase is in the melted meson phase, so the fluctuation 
spectrum is continuous, which means that the spectral "continuity" doesn't make 
much sense. However, in our A/" = 4 YM theory coupled to the M = 2 quark hy- 
permultiplet, it is known that the meson phase is continuously connected to a Higgs 
phase PT|. This marginal deformation does not cost any energy, and the baryon 
number density is kept to be zero. The instanton size modulus is a free parameter 
(that is, the squark VEV is a fiat direction of the theory). In this deformation, it 
was shown in Ref. |3^ that the discrete fluctuation spectrum is smoothly deformed. 
See Fig. 1 of Ref. |3^ . This phenomenon is analogous to the spectral quark-hadron 
continuity. 

It is well-known that color-flavor locking phase in QCD closely resembles the 
locking between spin and orbital symmetries found in the so-called "B-phase" of su- 
perfluid Helium 3, the setup we consider here therefore seems to be directly applicable 
in realizing this in string theory. One can study various thermodynamical properites 
and also consider topological defects e.g. vortices and study in such phase. Some 



interesting work relating D3D7 system with fermi-liquid can be found in Refs. [49 



Finally, it would be interesting to study a possible universality of the CFL at 
finite baryon density among holographic models. In the D4/D6 system considered in 
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Ref . p9[ , the dual field theory becomes effectively a pure bosonic Yang- Mills theory 
at low energy [50|. The phase structure of this system at finite temperature and 
baryon density was shown to have universal properties in Ref. Therefore, it is 
expected that when the baryon number density increases, the system becomes unsta- 
ble and some of the D4-branes would be pulled onto the D6 branes. In this case, the 
squarks condensation corresponds to an expansion of monopoles on the D6-branes, 
instead of the instantons. As mentioned in Sec. |3.1| , in the deconfinement phase, the 
baryon vertex is replaced by a flux while there is no probe brane description [^. On 
the other hand, in a confining phase, the baryon vertex does have a probe brane 
description, and the discussion in Sec. at zero temperature does not apply to the 
case. Therefore, it would be interesting to investigate the possibility of CFL in a 
confinement phase. 
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A. Check of Consistency for the Linearized Perturbation 



To complete the analysis of Sec. p.2| , we shall now check if this can be regarded as 
a small backreaction, so that our perturbative treatment for solving the equations 
of motion of the supergravity is guaranteed. The second term of (|3.13| ) suggests 



that the nonzero F3 will again backreact the F5 flux. We examine that this 

backreaction does not spoil the original flux configuration (|2.3|) too much. To this 
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end, we compare the second term of ( p.l3| ) with the kinetic term 

-1 



d''xy^o\F,\'. (A.l) 
We are only interested in order of magnitudes. Solving ( |3.17 ), we obtain 



Using this and (|3.1| ) (|2.3| ), we evaluate the second term of ( p.l3 



as 



/" F5 A ^2 A F3 ~ [ d^xdrf^dVt^ rfgla'^d^ . (A.3) 
4^10 J J 

On the other hand, the F5 kinetic term ([A.l| ) with the flux solution (|2.3| ) gives 

^ / t^'°Xv^|F5|2 ~ / d^xdr^dn, rlg-^a'-' . (A.4) 

Requiring ( |A.3|) being much smaller than (|A.4|) , we obtain 

gta'^'d^ < rl . (A.5) 

This means that, for the backreaction to the 5-form flux F5 to be small, we need to 
work in this region for rg. 

On the other hand, we made the assumption r <^ tq to simplify the source term 
to get (|3.12|) . Around the tip, we have a relation rg = + y(r)^ ~ r^(l + ?/'(0)), so 
this assupmtion translates to the condition rg(d^ — c^)/d^ ^ Tq which is equivalent 
to 

rl « «'8^2^7MV(d2 - c2)2 . (A.6) 

Therefore, in order to have a region for rg which satisfies the two requirements 
( [A75D and (^), we need 

^^.AT; « dV(d2 - c^) . (A.7) 

With ( |2.16| ) and ( |2.17| ), this condition is met if we are close to the critical chemical 
potential, 

/i — m-C/i,m. (A. 8) 

Throughout this paper, we are working in this regeme. 

Note that when d^ — c^ <^ d^ with which (|A.7|) is satisfied, the D7-brane spike be- 
comes very narrow, and the spike can be well-approximated by fundamental strings. 
This means that dilaton backreaction can be safely neglected. The backreaction to 
the metric is suppressed by 1/Nc and also neglected. 
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